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Finite temperature SU(3) gauge theory is studied on anisotropic lattices using the standard 
plaquette gauge action. The equation of state is calculated on 16 3 x 8, 20 3 x 10 and 24 3 x 12 lattices 
■ with the anisotropy £ = a 3 /a t = 2, where a a and a t are the spatial and temporal lattice spacings. 

Unlike the case of the isotropic lattice on which Nt = 4 data deviate significantly from the leading 
scaling behavior, the pressure and energy density on an anisotropic lattice are found to satisfy well 
the leading l/JVf scaling from our coarsest lattice, Nt/£ = 4. With three data points at Nt/£ = 4, 
5 and 6, we perform a well controlled continuum extrapolation of the equation of state. Our results 
in the continuum limit agree with a previous result from isotropic lattices using the same action, 
but have smaller and more reliable errors. 
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I. INTRODUCTION 



Study of lattice QCD at finite temperatures is an important step towards clarification of the dynamics of the quark 
gluon plasma which is believed to have formed in the early Universe and is expected to be created in high energy heavy 
ion collisions (!]]. In order to extract predictions for the real world from results obtained on finite lattices, we have 
to extrapolate lattice data to the continuum limit of vanishing lattice spacings. Because of the large computational 
JL , demands for full QCD simulations, continuum extrapolations of thermodynamic quantities have so far been attempted 
only in SU(3) gauge theory, i.e., in the quenched approximation of QCD, where the influence of dynamical quarks 
is neglected. Two studies using the standard plaquette gauge action ||] and a renormalization-group (RG) improved 
L* ■ gauge action || have found the pressure and energy density consistent with each other in the continuum limit. 
. £^ | In full QCD with two flavors of dynamical quarks, thermodynamic quantities on coarse lattices have been found 
to show large lattice spacing dependence @-|(J. For a reliable extrapolation to the continuum limit, data on finer 
lattices are required. With conventional isotropic lattices, this means an increase of the spatial lattice size to keep 
the physical volume close to the thermodynamic limit. Full QCD simulations on large lattices are still difficult with 
the current computer power. A more efficient method of calculation is desirable. Even in the quenched case, we note 
that continuum extrapolations of equation of state have been made using only two lattice spacings ^ n order 

to reliably estimate systematic errors from the extrapolations, more data points are needed. Therefore, an efficient 
method is called for also in quenched QCD. 

Recently, anisotropic lattices have been employed to study transport coefficients and temporal correlation functions 
in finite temperature QCD |t] 6|. In these studies, anisotropy was introduced to obtain more data points for temporal 
correlation functions. 

In this paper, we show that anisotropic lattices provide also an efficient calculation method for thermodynamic 
quantities. The idea is as follows. Inspecting the free energy density of SU(3) gauge theory in the high temperature 
Stephan-Boltzmann limit, the leading discretization error from the temporal direction is found to be much larger 
than that from each of the spatial directions. Hence, choosing £ = a s /a t larger than one, where a s and at are the 
spatial and temporal lattice spacings, cutoff errors in thermodynamic quantities will be efficiently reduced without 
much increase in the computational cost. From a study of free energy density in the high temperature limit, we find 
that £ = 2 is an optimal choice for SU(3) gauge theory. This improvement also makes it computationally easier to 
accumulate data for more values of temporal lattice sizes for the continuum extrapolation. 

As a first test of the method, we study the equation of state (EOS) in SU(3) gauge theory. On isotropic lattices, 
discretization errors in the EOS for the plaquette action are quite large at the temporal lattice size Nt — 4. The data 
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at this value of N t deviate significantly from the leading 1/iV 2 scaling behavior, F(T)\ N = ^(T)| continuum + c f/N^, 
where F is a thermodynamic quantity at a fixed temperature T. So far, continuum extrapolations of the EOS have 
been made using results at Nt — 6 and 8. On anisotropic lattices with £ = 2, we find that the discretization errors 
in the pressure and energy density are much reduced relative to those from isotropic lattices with the same spatial 
lattice spacing. Furthermore, we find that the EOS at Nt/£ = 4, 5 and 6 follow the leading scaling behavior 

remarkably well. Therefore, a continuum extrapolation can be reliably carried out. Since the total computational 
cost is still lower than that for an Nt — 8 isotropic simulation, we can achieve a higher statistics as well, resulting in 
smaller final errors. 

In Sec. |n|, we study the high temperature limit of SU(3) gauge theory on anisotropic lattices to see how £ appears 
in the leading discretization error for the EOS. From this study, we find that £ = 2 is an optimum choice for our 
purpose. We then perform a series of simulations on £ = 2 anisotropic lattices. Our lattice action and simulation 



parameters arc described in Sec. III. Sec. IV is devoted to a calculation of the lattice scale through the string tension. 
The critical temperature is determined in Sec. [y|. Our main results are presented in Sees. VI and VII, where the 



pressure and ener gy density are calculated and their continuum extrapolations are carried out. A brief summary is 



given in Sec. VIII 




II. HIGH TEMPERATURE LIMIT 

In the high temperature limit, the gauge coupling vanishes due to asymptotic freedom, and SU(3) gauge theory 
turns into a free bosonic gas. In the integral method JlCfl which we apply in this study, the pressure p is related to 
the free energy density / by p = — f for large homogeneous systems. Therefore, in the high temperature limit, the 
energy density e is given by e = 3p = —3/. The value of / in the high temperature limit has been calculated in ]TT] , p~2[ . 
Normalizing e by the Stephan-Boltzmann value in the continuum limit, we find 

e i 5 + 3^ 2 / tt \ 2 91 + 210£ 2 + 99£ 4 / it \ 4 ~ /V tt \ 8 \ 
- B = 1 + —{W t ) + 1680 U) +0 (U) j (1) 

for spatially large lattices. Substituting £ = 1 in Eq. ([!]), we recover the previous results for isotropic lattices [ pi) . 
When we alternatively adopt the derivative method (operator method) jll| to define the energy density, we obtain 

(2) 

In both formulae, the leading discretization error is proportional to l/A" 2 . 

In the leading 1/iV 2 term of Eq. ([!]) (or Eq. (0)), the term proportional to £ 2 represents the discretization error 
from finite lattice spacings a s in the three spatial directions. We find that the temporal cutoff at leads to 5/8 (or 1/2) 
of the leading discretization error at £ = 1, while the spatial cutoff a s contributes only 1/8 (or 1/6) from each of the 
three spatial directions. 

Since a reduction of the lattice spacing in each direction separately causes an increase of the computational cost 
by a similar magnitude, a reduction of a t is much more efficient than that of a s in suppressing lattice artifacts in 
thermodynamic quantities. Making the anisotropy £ = a s /at too large is, however, again inefficient because the 
spatial discretization errors remain even in the limit of £ = oo. A rough estimate for the optimum value of £ is 
given by equating the discretization errors from spatial and temporal directions, £ = \fh ~ 2.24 from Eq. ([!]), and 
£ = V3 w 1.73 from Eq. (||). More elaborate estimations considering the balance between the computational cost as 
a function of the lattice size and the magnitude of discretization errors including higher orders of 1 /N t lead to similar 
values of £. 

Based on these considerations, we adopt £ = 2 for simulations of SU(3) gauge theory in the present work. An 
even number for £ is attractive also for the vectorization/parallelization of the simulation code which is based on an 
even-odd algorithm, since we can study the case of odd AT t /£ without modifying the program. 

III. DETAILS OF SIMULATIONS 

A. Action 

We employ the plaquette gauge action for SU(3) gauge theory given by 
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S G [U]=pQ-Q s + toQt), (3) 
where £o is the bare anisotropy, f3 = 6/g 2 with go the bare gauge coupling constant, and 

q s = a - , Qt = E c 1 - p ^ - ( 4 ) 

n,(tj) n,i 

with P[j,v(n) = iRe Tr U^ v (n) the plaquette in the (//, ^) plane at site n. Anisotropy is introduced by choosing £ 7^ 1. 

Due to quantum fluctuations, the actual anisotropy £ = a s /a t deviates from the bare value £o- We define the 
renormalization factor r)(/3,£) for £ by 

The values of rj(@,£) can be determined non-perturbatively by matching Wilson loops in temporal and spatial direc- 
tions on anisotropic lattices fl^-|l6f|. For our simulation, we calculate £o(/?>£ = 2) using r]((3,£) obtained by Klassen 
for the range 1 < £ < 6 and 5.5 < < 00 Eg]: 



1\ 1 + Qigg 

6 l + a 5§'' 



^0 = l+(l-7)^^§5o 2 , (6) 



where a = -0.77810, a x = -0.55055 and 



„ _ 1.002503£ 3 + 0.39100£ 2 + 1.47130£ - 0.19231 

Vl ^' ~ £ 3 + 0.26287£ 2 + 1.59008£ - 0.18224 ' ^ ^ 



B. Simulation parameters 

The main runs of our simulations are carried out on £ = 2 anisotropic lattices with size x N t = 16 3 x 8, 20 3 x 10 
and 24 3 x 12. For N t — 8, we make additional runs on 12 3 x 8 and 24 3 x 8 lattices to examine finite size effects. The 
zero-temperature runs are made on iV 3 x £N S lattices with £ = 2. The simulation parameters of these runs which 
cover the range T/T c ~ 0.9-5.0 are listed in Table |[ To determine precise values for the critical coupling, longer runs 
around the critical points are made at the parameters compiled in Table |Tl[ 

For the main runs, the aspect ratio L S T = (N s a s )/ (N t at) is fixed to 4, where L s = N s a s is the spatial lattice size in 



physical units. This choice is based on a study of finite spatial volume effects presented in Sec. VI, where it is shown 
that, for the precision and the range of T/T c we study, finite spatial volume effects in the EOS are sufficiently small 
with L S T > 4. 

Gauge configurations are generated by a 5-hit pseudo heat bath update followed by four over-relaxation sweeps, 



which we call an iteration. As discussed in Sec. |VI], the total number of iterations should be approximately proportional 
to 7V t 6 to keep an accuracy for the EOS. After thermalization, we perform 20,000 to 100,000 iterations on finite- 
temperature lattices and 5,000 to 25,000 iterations on zero-temperature lattices, as compiled in Table |[ At every 
iteration, we measure the spatial and temporal plaquettes, P ss and P st . Near the critical temperature, we also measure 
the Polyakov loop. The errors are estimated by a jack-knife method. The bin size for the jack-knife errors, listed in 
Table |, is determined from a study of bin size dependence as illustrated in Fig. [I]. The results for the plaquettes are 
summarized in Tables III-^. 



IV. SCALE 
A. Static quark potential 

We determine the physical scale of our lattices from the string tension, which is calculated from the static 
quark-antiquark potential at zero temperature. To calculate the static quark potential, we perform additional zero- 



temperature simulations listed in Table VI. The static quark potential V(R) is defined through 

W(R,f) =C{R)e- v{k)f/ t, (8) 
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where W(R, f) is the Wilson loop in a spatial-temporal plane with the size Ra s x Tat- We measure Wilson loops at 
every 25 iterations after thermalization. In order to enhance the ground state signal in (||), we smear the spatial links 
of the Wilson loop Jl7],|l8|]. Details of the smearing method are the same as in Ref. |Q. We determine the optimum 
smearing step N opt which maximizes the overlap function C(R) under the condition C(R) < 1. Following Ref. [ Ell , 
we study a local effective potential defined by 



V eff (R,T)=tlog 



( W(R,f) 
\W(R,f +1) 



(9) 



which tends to V(R) at sufficiently large T. The reason to adopt Eq. (^) instead of the fit result from Eq. (||) is to 



perform a correlated error analysis directly for the potential parameters. The optimum value of T, listed in Table VII 
is obtained by inspecting the plateau of V e ff(R, T) at each (3. 

We perform a correlated fit of V(R) = V e jj(R,T op t) with the ansatz p3), 



V(R) = V a +<jR-el- + l(lr- 
R \R 



(10) 



Here, 



is the lattice Coulomb term from one gluon exchange 

d 3 k cos(k • R) 



" 1 " 






= Air J 


.R. 





(11) 



which is introduced to approximately remove terms violating rotational invariance at short distances. The coefficient 
I is treated as a free parameter. 

The fit range [R m in, Rmax] for R is determined by consulting the stability of the fit. Our choices for R m in are given 
in Table VII. We confirm that the fits and the values of the string tension are stable under a variation of R m in- The 
string tension is almost insensitive to a wide variation of R max . Hence R m ax is chosen as large as possible so far as 
the fit is stable and the signal is not lost in the noise. With this choice for the fit range, we obtain fit curves which 
reproduce the data well. 



Our results for the potential parameters are summarized in Table VH. The error includes the jack-knife error with 
bin size one (25 iterations) and the systematic error from the choice of R m in estimated through a difference under the 
change of R m i n by one. We confirm that increasing the bin size to two gives consistent results on 16 3 x 32 lattices, 
while, on 24 3 x 48 lattices, correlated fits with bin size two become unstable due to insufficient number of jackknife 
ensembles. 



B. String tension 

We interpolate the string tension data using an ansatz proposed by Allton [pl| , 



1 + c 2 a(/3) 2 + c 4 a(/3) 4 
co 



where /(/3) is the two-loop scaling function of SU(3) gauge theory, 

f(P) = Hr ex P[ 



13 J 126 J ' 
, 11 102 

60 -16^' 6l= (lW ( } 

andaC8)=/09V/(/3 = 6.O). 

From Table VII, we find that the values for a s ^fa are insensitive to the spatial lattice volume to the present 
precision. Using data marked by star (*) in Table VII, we obtain the best fit at 

c = 0.01171(41), c 2 = 0.285(79), c 4 = 0.033(30), (14) 

with x 2 /Ndf — 1.77. The string tension data and the resulting fit curve are shown in Fig. together with those 
from isotropic lattices |22] . 
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V. CRITICAL TEMPERATURE 



We define the critical gauge coupling (3 c (N t , N s ) from the location of the peak of the susceptibility Xrot for a Z(3)- 
rotated Polyakov loop. The simulation parameters for the study of j3 c are compiled in Table O. Th e /^-d ependence of 



Xrot is calculated using the spectral density method |23j . The results for j3 c are compiled in Table |VIII . 

To estimate the critical temperature, we have to extrapolate P c (N t , N s ) to the thermodynamic limit and to the 
continuum limit. We perform the extrapolation to the thermodynamic limit using a finite-size scaling ansatz, 

p c (N t ,N s )=p c (N t ,<x>)-h(J^) . (15) 

for first order phase transitions. From the data for (3 C on anisotropic 12 3 x 8, 16 3 x 8 and 24 3 x 8 lattices with £ = 2, 
we find h — 0.031(16) for N t /£, = 4, as shown in Fi g, pj . In a previous study on isotropic lattices, h was found to be 
approximately independent of N t for N t = 4 and 6 |24J . We extract (3 c {N t , oo) adopting h = 0.031(16) for all N t . 
The critical temperature in units of the string tension is given by 

7^ = N t (a ay /F) (/J c (iV t ,oo)) (16) 



using the fit result for Eq. (|l2|). The values of T c j ^[5 are summarized in Fig. |Z] and Table VIII . The dominant part 
of the errors in T c j 'yfa is from the Allton fit for the string tension. 

Finally we extrapolate the results to the continuum limit assuming the leading l/iV t 2 scaling ansatz, 

F\Nt = ^ I continuum + jy2 (^) 

with F = T c j \fo. The extrapolation is shown in Fig. |J. In the continuum limit, we obtain 

= 0.635(10) (18) 

V' 7 

from the £ = 2 plaquette action. 

In Fi g. El we also plot the results obtained on isotropic lattices using the plaquette action |^5j and the RG-improved 
action pOH. Our value of T c j^fa in the continuum limit is consistent with these results within the error of about 
2%. A more precise comparison would require the generation and analyses of potential data in a completely parallel 
manner, because, as discussed in ||, numerical values of T c /y/a at a few percent level sensitively depend on the 
method used to determine the string tension. We leave this issue for future studies. 

VI. PRESSURE 
A. Integral method 

We use the integral method to calculate the pressure [0. This method is based on the relation p = — f = 
(T/V) log Z(T, V) satisfied for a large homogeneous system, where V — £ 3 is the spatial volume of the system in 
physical units and Z is the partition function. Rewriting logZ = f d{3j?jj^, the pressure is given by 



P_ 



= fdp&Stf), (19) 

00 J/3 



with 



iVA 4 1 dlogZ 



£ J NfN t dp 

For our anisotropic gauge action (^J), the derivative of log Z is given by 



(20) 



aiogz _/s G \ , pdtoifi Q / _ igh, . _ = Q contribution) (21) 



d/3 \ & I " dp V 
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We use symmetric iV 3 x £N S lattices to calculate the T = contribution. For a sufficiently small fto, p(Po) can be 
neglected. 

In order to keep the same accuracy of AS" for the same physical lattice volume L 3 in units of the temperature 
T, the statistics of simulations should increase in proportion to (£,(N t / 1;) 4 ) 2 / (N^ N t ) cx N 4 /^ 3 . Here, the first factor 
arises from £(-?Vt/£) 4 in Eq. (pp|), and the second factor l/(N 3 N t ) from a suppression of fluctuations due to averaging 
over the lattice volume. Taking into account the autocorrelation time which is proportional to N 2 , the number of 
iterations should increase as ~ N®. 

Integrating AS* in f3 using a cubic spline interpolation, we obtain the pressure. For the horizontal axis, we use the 
temperature in units of the critical temperature, 

T c (<W?)(/3)' 

The errors from numerical integration are estimated by a jack-knife method in the following way ||. Since simulations 
at different (3 are statistically independent, we sum up all the contributions from fit smaller than (3 corresponding to 
the temperature T by the naive error-propagation rule, Sp(T) = 6pi(T) 2 , where Sip(T) at each simulation point 

Pi is estimated by the jack-knife method. 



B. Finite spatial volume effects 

We first study the effects of finite spatial volume on the EOS. In Fig. [|, we show the results for AS at iV t /£ = 8/2 
with the aspect ratio L S T = N s £/N t — 3, 4 and 6 which correspond to N s = 12, 16 and 24, respectively. Integrating 
AS in /3, we obtain Fig. ||for the pressure. We find that the data at L S T = 3 is affected by sizable finite volume effects 
both at T ~ T c and at high temperatures. On the other hand, for the range of T/T c we study, the pressure does not 
change when the aspect ratio is increased from L S T — 4 to 6, indicating that the conventional choice L S T — 4 is safe 
with the present precision of data. Hence, we choose L S T — 4 for our studies of lattice spacing dependence. Results 
for AS at L S T = 4 with various N t are given in Fig. |?]. Integrating the data using a cubic spline interpolation, as 
shown in the figures, we obtain the pressure plotted in Fig. ||. 



C. Continuum extrapolation 

We now extrapolate the pressure to the continuum limit using the leading order ansatz of Eq. (|l7|) . Figure ^| shows 
the pressure at T/T c = 1.5, 2.5 and 3.5 as a function of (£/iV t ) 2 (filled circles). For comparison, results from isotropic 
lattices using the plaquette action || (open circles) and the RG-improved action || (open squares) are also plotted. 
For the £ = 1 plaquette data, we adopt the results of a reanalysis made in Ref. || to commonly apply the scale from 
the Allton fit of the string tension and also the same error estimation method. 

The advantage of using anisotropic lattices is apparent from Fig. Q On the coarsest lattice Nt/£ = 4, finite lattice 
spacing errors at £ = 2 are much smaller than those at £ = 1 with the same plaquette action. The pressure at 
T = 2.5T C , for example, on the isotropic 16 3 x 4 lattice is larger than its continuum limit by about 20%, while the 
deviation is only 5% on the corresponding 16 3 x 8 lattice with £ = 2. Furthermore, with the anisotropic £ = 2 data, 
the leading l/iV t 2 term describes the data well even at Nt/£ = 4 (the right-most point). Therefore, we can confidently 
perform an extrapolation to the continuum limit using three data points. In the case of the isotropic plaquette action, 
in contrast, the continuum extrapolation had to be made with only two data points at N t /£ — 6 and 8. In the 
continuum limit, our results for £ = 2 are slightly smaller than those from the isotropic plaquette action, but the 
results are consistent with each other within the error of about 5% for the results from the isotropic action. It is 
worth observing that the £ = 2 results have smaller and more reliable errors of 2-3%. 

In order to quantitatively evaluate the benefit of anisotropic lattices, we compare the computational cost to achieve 
comparable systematic and statistical errors on isotropic and £ = 2 anisotropic lattices. Choosing T — 2.5T C as 
a typical example, we find that the deviation of the pressure from the continuum limit (i.e., the magnitude of 
the systematic error due to finite lattice cutoffs) is comparable between the isotropic 32 3 x 8 and our £ = 2 
anisotropic 20 3 x 10 lattices, i.e., p/T 4 = 1.390(26) on a 32 3 x 8 lattice and p/T 4 = 1.381(13) on a 20 3 x 10 
lattice, both lattices having the same spatial size N s a s = 1.6/T C . The number of configurations to achieve these 
statistical errors are 20,000-40,000 iterations for £ = 1 and 50,000 for £ = 2, respectively. Therefore, for the same 
statistical error, the relative computational cost for a £ = 2 lattice over that for £ = 1 is conservatively estimated 
as ((20 3 x 10) x 50000) / ((32 3 x 8) x 4 x 20000) w 1/5, showing a factor 5 gain in the computational cost for the 
anisotropic calculation in this example. 
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In Fig. ||we also note that the results from the RG- improved action on isotropic lattices are higher by 7-10% (about 
2a) than those from the present work in the continuum limit. A possible origin of this discrepancy is the use of the 
iV t /£ = 4 data of the RG- improved action, which show a large (about 20%) deviation from the continuum value. For 
a detailed test of consistency, we need more data points, say at -/V t /£ = 6, from the RG-improved action. 

Repeating the continuum extrapolation at other values of T/T c , we obtain Fig. [w| Our results show a quite slow 
approach to the high temperature Stcphan-Boltzmann limit, as reported also in previous studies on isotropic lattices 

§|- 



VII. ENERGY DENSITY 



by 



We calculate the energy density e by combining the results of p/T 4 with those for the interaction measure defined 



t — 'ip 



03 



T 4 ~ a& da s 



AS. 



(23) 



The QCD beta function on anisotropic lattice 



813 



do 



12&o 



is determined through the string tension a studied in Sec. IV B, 

(24) 



1 + c 2 a 2 + c 4 a 4 



£ 6 (6i/6 ) [3- 1 - 1 1 + 3c 2 a 2 + 5c 4 a 4 ' 



where the coefficients Ci are given in Eq. (|l2]). The error of the energy density is calculated by quadrature from the 
error of 3p and that for e — 3p, the latter being proportional to the error of AS. 

The results for the energy density are shown in Figs, [ll] and [l2|. As in the case of the pressure the leading scaling 
behavior is well followed by our £ = 2 data from N t /£, = 4, which allows us to extrapolate to the continuum limit 
reliably. The results for the energy density in the continuum limit are compared with the previous results in Fig. 
Our £ = 2 plaquette action leads to an energy density which is slightly smaller than, but consistent with that from 
the £ = 1 plaquette action, but is about 7-10% (about 2a) smaller than that from the £ = 1 RG action. More work 
is required to clarify the origin of the small discrepancy with the RG action. 



VIII. CONCLUSION 



We have studied the continuum limit of the equation of state in SU(3) gauge theory on anisotropic lattices with 
the anisotropy £ = a s /at — 2, using the standard plaquette gauge action. Anisotropic lattices are shown to be 
more efficient in calculating thermodynamic quantities than isotropic lattices. We found that the cutoff errors in the 
pressure and energy density arc much smaller than corresponding isotropic lattice results at small values of N t /t;. The 
computational cost for £ = 2 lattices is about 1/5 of that for £ = 1 lattices. We also found that the leading scaling 
behavior is well satisfied already from iV t /£ = 4, which enabled us to perform continuum extrapolations with three 
data points at N t /£, = 4, 5 and 6. The equation of state in the continuum limit agrees with that obtained on isotropic 
lattice using the same action, but have much smaller and better controlled errors. The benefit of anisotropic lattice 
demonstrated here will be indispensable for extraction of continuum predictions for the equation of state, when we 
include dynamical quarks. 
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lattice P bin size # of iter. 

12 3 x 8 5.73-6.80 1600 40 000 

16 3 x 8* 5.74-6.80 800 20 000 

24 3 x 8 5.75-6.80 400 10 000 

20 3 x 10* 5.86-6.98 2000 50 000 

24 3 x 12* 5.95-7.20 4000 100 000 

12 3 x 24 5.74-6.80 400 10 000 

16 3 x 32* 5.74-6.80 200 5 000 

20 3 x 40* 5.86-6.98 500 12 500 

24 3 x 48 5.75-5.90 100 2 500 

24 3 x 48* 5.95-7.20 1000 25 000 

TABLE I. Simulation parameters. Main runs are marked by star (*). 



lattice /3 bin size # of iter. 

12 3 x 8 5.790, 5.791 8000 80 000 

16 3 x 8 5.790 , 5.792 4000 40 000 

24 3 x 8 5.791, 5.792 4000 40 000 

20 3 x 10 5.903 , 5.907 5000 50 000 

24 3 x 12 6.004, 6.006 10000 100 000 



TABLE II. Simulation parameters for determination of critical couplings. 
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16 3 x 8 16 3 x 32 



p 


Co 


p 

J ss 


Pst 


p 


Pst 


5.740 


1.66279318 


0.448467(31) 


0.679985(12) 


0.448490(28) 


0.679979(11) 


5.750 


1.66473308 


0.450693(24) 


0.681412(11) 


0.450641(21) 


0.681384(8) 


5.760 


1.66664410 


0.452784(33) 


0.682783(13) 


0.452731(22) 


0.682747(9) 


5.770 


1.66852693 


0.454935(29) 


0.684175(13) 


0.454758(24) 


0.684090(9) 


5.780 


1.67038223 


0.457024(53) 


0.685533(22) 


0.456720(21) 


0.685392(8) 


5.788 


1.67184708 


0.459186(116) 


0.686823(49) 


0.458272(30) 


0.686419(11) 


5.790 


1.67221065 


0.459930(109) 


0.687240(48) 


0.458678(26) 


0.686679(11) 


5.792 


1.67257316 


0.460517(104) 


0.687578(45) 


0.459056(22) 


0.686929(9) 


5.800 


1.67401280 


0.462698(75) 


0.688873(33) 


0.460586(22) 


0.687949(9) 


5.805 


1.67490422 


0.463825(34) 


0.689587(15) 


0.461565(21) 


0.688588(9) 


5.810 


1.67578929 


0.464912(40) 


0.690278(17) 


0.462446(20) 


0.689181(9) 


5.820 


1.67754071 


0.466746(21) 


0.691520(10) 


0.464241(17) 


0.690383(6) 


5.830 


1.67926762 


0.468486(24) 


0.692704(10) 


0.466022(21) 


0.691578(9) 


5.840 


1.68097058 


0.470122(18) 


0.693839(8) 


0.467707(24) 


0.692722(9) 


5.880 


1.68755324 


0.476195(15) 


0.698142(7) 


0.474205(17) 


0.697145(7) 


5.900 


1.69071395 


0.478994(18) 


0.700156(9) 


0.477282(22) 


0.699255(9) 


5.950 


1.69826359 


0.485606(15) 


0.704933(7) 


0.484390(18) 


0.704199(7) 


6.000 


1.70535029 


0.491774(15) 


0.709406(6) 


0.490955(20) 


0.708801(9) 


6.100 


1.71830738 


0.503237(14) 


0.717652(6) 


0.502986(14) 


0.717230(5) 


6.200 


1.72987892 


0.513833(11) 


0.725175(6) 


0.513839(14) 


0.724837(5) 


6.300 


1.74029271 


0.523743(10) 


0.732106(4) 


0.523915(15) 


0.731827(7) 


6.400 


1.74972820 


0.533075(11) 


0.738552(4) 


0.533401(9) 


0.738316(3) 


6.500 


1.75832876 


0.541970(13) 


0.744586(5) 


0.542362(8) 


0.744378(5) 


6.600 


1.76621035 


0.550391(8) 


0.750250(3) 


0.550854(10) 


0.750058(4) 


6.700 


1.77346785 


0.558485(9) 


0.755608(4) 


0.558959(9) 


0.755427(4) 


6.800 


1.78017964 


0.566215(12) 


0.760672(4) 


0.566716(8) 


0.760501(4) 




TABLE III. 


Plaquette expectation values 


on 16 3 x 8 and 16 3 


x 32 lattices with £ = 2. 





10 



20 3 x 10 20 3 x 40 



p 


Co 


p 

J ss 


Pst 


p 


Pst 


5.86288916 


1.68478116 


0.4715286(90) 


0.6953072(38) 


0.4715194(98) 


0.6953039(38) 


5.87 


1.68594094 


0.4726803(97) 


0.6960907(37) 


0.4726453(79) 


0.6960771(33) 


5.88583578 


1.68848420 


0.4752043(113) 


0.6978062(52) 


0.4751072(93) 


0.6977655(41) 


5.90 


1.69071395 


0.4775533(342) 


0.6993698(144) 


0.4772612(79) 


0.6992430(33) 


5.91 


1.69226327 


0.4793349(340) 


0.7005240(144) 


0.4787235(65) 


0.7002573(30) 


5.92 


1.69379248 


0.4809915(113) 


0.7016191(50) 


0.4801832(57) 


0.7012665(26) 


5.93084722 


1.69542899 


0.4826008(89) 


0.7027227(39) 


0.4817182(78) 


0.7023359(35) 


5.94 


1.69679224 


0.4838962(61) 


0.7036250(30) 


0.4830113(60) 


0.7032314(30) 


5.96 


1.69971645 


0.4865820(62) 


0.7055225(30) 


0.4857427(62) 


0.7051382(32) 


5.98 


1.70256818 


0.4891795(54) 


0.7073650(25) 


0.4883883(83) 


0.7069900(34) 


5.9961937 


1.70482605 


0.4912217(55) 


0.7088160(30) 


0.4904832(71) 


0.7084591(30) 


6.0793640 


1.71575557 


0.5010417(44) 


0.7158270(31) 


0.5005840(62) 


0.7155576(27) 


6.17716193 


1.72734556 


0.5116532(54) 


0.7233550(25) 


0.5114357(43) 


0.7231598(22) 


6.28582916 


1.73888020 


0.5225991(56) 


0.7310157(21) 


0.5225280(53) 


0.7308687(21) 


6.40118969 


1.74983517 


0.5334631(32) 


0.7385009(19) 


0.5334926(43) 


0.7383839(17) 


6 51881026 


1 75986308 


0.5438681(48) 


0.7455581(19) 


5439702(40) 


7454657(19) 


6.63417079 


1.76875624 


0.5535144(38) 


0.7520032(19) 


0.5536476(51) 


0.7519204(23) 


6.74283803 


1.77640579 


0.5621461(45) 


0.7576970(23) 


0.5623098(36) 


0.7576251(14) 


6.84063596 


1.78276647 


0.5695876(32) 


0.7625475(17) 


0.5697626(34) 


0.7624799(11) 


6.92380626 


1.78783002 


0.5756793(33) 


0.7664882(18) 


0.5758587(31) 


0.7664206(16) 


6.98915275 


1.79160648 


0.5803248(35) 


0.7694702(14) 


0.5805094(41) 


0.7694057(18) 




TABLE IV. 


Plaquette expectation values on 20 3 x 10 and 20 3 


x 40 lattices with £ = 2. 





24 3 x 12 24 3 x 48 



p 


Co 


p 

J ss 


Pst 


p 

1 ss 


Pst 


5.95 


1.69826359 


0.4843851(27) 


0.7041916(13) 


0.4843789(45) 


0.7041883(19) 


5.98 


1.70256818 


0.4884099(39) 


0.7070003(19) 


0.4883825(35) 


0.7069880(15) 


6.00 


1.70535029 


0.4911005(118) 


0.7088537(50) 


0.4909663(38) 


0.7087977(14) 


6.01 


1.70671610 


0.4924924(104) 


0.7097962(43) 


0.4922291(37) 


0.7096838(15) 


6.02 


1.70806552 


0.4938053(64) 


0.7107011(32) 


0.4934718(30) 


0.7105575(13) 


6.03 


1.70939887 


0.4950807(40) 


0.7115881(16) 


0.4947043(36) 


0.7114232(17) 


6.04 


1.71071646 


0.4963132(30) 


0.7124510(16) 


0.4959199(32) 


0.7122791(13) 


6.07 


1.71457763 


0.4998634(27) 


0.7149595(10) 


0.4994891(31) 


0.7147889(15) 


6.08 


1.71583512 


0.5010194(19) 


0.7157747(6) 


0.5006575(31) 


0.7156082(13) 


6.10 


1.71830738 


0.5032879(22) 


0.7173807(10) 


0.5029551(29) 


0.7172208(13) 


6.15 


1.72425080 


0.5087787(26) 


0.7212576(10) 


0.5085106(19) 


0.7211154(12) 


6.20 


1.72987892 


0.5140368(26) 


0.7249549(12) 


0.5138372(20) 


0.7248368(8) 


6.30 


1.74029271 


0.5240287(21) 


0.7319188(8) 


0.5239220(23) 


0.7318284(10) 


6.40 


1.74972820 


0.5334259(25) 


0.7383798(11) 


0.5333873(23) 


0.7383125(9) 


6.60 


1.76621035 


0.5508062(15) 


0.7501014(7) 


0.5508372(22) 


0.7500563(9) 


6.80 


1.78017964 


0.5666348(15) 


0.7605281(6) 


0.5667010(21) 


0.7604924(9) 


7.00 


1.79221720 


0.5811933(20) 


0.7699251(8) 


0.5812721(12) 


0.7698933(6) 


7.20 


1.80273290 


0.5946688(17) 


0.7784726(9) 


0.5947568(18) 


0.7784435(8) 




TABLE V. 


Plaquette expectation values 


on 24 3 x 12 and 24 3 


x 48 lattices with i = 2. 
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p 


lattice 


AT 

JVopt 


or com. 


5.7 


16 3 x 32 


3 


800 


5.8 


16 3 x 32 


5 


800 


5.9 


16 3 x 32 


6 


800 


6.0 


16 3 x 32 


8 


600 




24 3 X 48 


8 


100 


6.1 


16 3 x 32 


10 


400 


6.3 


16 3 x 32 


16 


300 




24 3 X 48 


20 


100 


6.5 


24 3 X 48 


30 


100 



TABLE VI. Simulation parameters for static quark potential at zero temperature. 



p 


lattice 


a s \fcr 


L s [tm] 


f 


Rinin 


V 


e 


I 


X 2 /Ndf 


5.7 


16 3 x 32* 


0.4794(66) 


3.49 


5 


V5 


0.677(36) 


0.305(50) 


0.934(122) 


5.81 


5.8 


16 3 x 32* 


0.3804(24) 


2.77 


6 


VE 


0.720(11) 


0.326(16) 


0.647(49) 


3.07 


5.9 


16 3 x 32* 


0.3190(18) 


2.32 


7 




0.688(7) 


0.284(11) 


0.501(43) 


3.20 


6.0 


16 3 x 32 


0.2667(21) 


1.94 


8 


Vg 


0.685(8) 


0.283(14) 


0.396(73) 


0.93 




24 3 X 48* 


0.2611(31) 


2.85 


8 


v/6 


0.699(11) 


0.310(19) 


0.565(82) 


2.05 


6.1 


16 3 x 32* 


0.2224(20) 


1.61 


8 


2^2 


0.686(6) 


0.297(13) 


0.375(61) 


1.97 


6.3 


16 3 x 32 


0.1656(19) 


1.20 


9 




0.653(5) 


0.281(9) 


0.239(67) 


0.95 




24 3 X 48* 


0.1661(20) 


1.81 


9 


v& 


0.657(5) 


0.294(9) 


0.323(68) 


1.72 


6.5 


24 3 X 48* 


0.1242(21) 


1.35 


9 


Vg 


0.622(3) 


0.279(6) 


0.247(47) 


1.75 



TABLE VII. Results for the potential parameters on £ = 2 anisotropic lattices with the plaquette action. The spatial lattice 
size L s is computed using ^fa — 440 MeV. 



iV 3 x N t 


12 3 x 8 16 3 x 8 


24 3 x 8 


20 3 x 10 


24 3 x 12 


Pc{N t ,N s ) 


5.79037(40) 5.79081(54) 


5.79138(31) 


5.90494(92) 


6.00464(67) 


p c (N t ,oo) 


5.79149(34) 




5.90543(116) 


6.00512(91) 




0.6402(39) 




0.6392(39) 


0.6364(75) 


TABLE VIII. 


Critical coupling and temperature on 


anisotropic d; = 2 lattices. 


Results for T c j\fa are 


obtained in the 



thermodynamic limit. 



12 



0.00005 




0.00000 




1000 1500 

Bin size 

FIG. 1. Typical bin size dependence of jack-knife errors for AS. 
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FIG. 3. Finite-size scaling of f3 c for Nt/£ = 4 on £ = 2 anisotropic lattices. 
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FIG. 5. Spatial lattice volume dependence in AS at Nt/£ = 4 on N s = 12, 16 and 24 lattices with £ = 2. 




2.5 3 

FIG. 6. Spatial volume dependence of the pressure p/T 4 on £ = 2 anisotropic lattices with N t /£ = 4. 



15 




CO 



5.8 


6 


6.2 6.4 6.6 6.8 

P 


? \ 

| X 




Nt= 10 


j 


• 




/ 









5.8 6 6.2 6.4 6.6 6.8 

P 



16 



i \ 

\ 








Nt= 12 


i \ 

* 






























4 














6 6.2 6.4 6.6 6.8 7 7.2 

FIG. 7. AS on JV t /| = 4, 5 and 6 lattices with £ = 2. 




FIG. 8. Pressure p/T 4 on £ = 2 anisotropic lattices. 
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FIG. 10. Pressure p/T 4 in the continuum limit. 
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FIG. 11. e/T 4 on anisotropic 16 3 x 8, 20 3 x 10 and 24 3 x 12 lattices with £ = 2. 
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